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Scalar color octets are generic signals for new physics at LHC energies. We examine tiieir pair 
production at the LHC to next-to-leading order QCD. This computation serves as another test of the 
fully automized MadGolem framework. We find large NLO production rates and sizeable quantum 
effects which depend on the sgluon mass. The shift in the sgluon distributions is mild and in good 
agreement with a multi-jet merging calculation. 

I. INTRODUCTION 



o 

^Nl Sgluons [H [2] are a type of scalar color octet states which arises in a variety of extensions of the Standard Model [3] . 

; I They can be fundamental or composite degrees of freedom. In extended supersymmetric models like the i?-symmetric 

Ci MSSM [llllllS] or A/" = l/Af — 2 hybrid models |2], sgluons emerge as scalar partners of a Dirac gluino. More generally, 

^> sgluons appear to be ubiquitous in models of supersymmetry breaking [3] . Compositeness models include fermion fields 
which transform under a confining gauge group — some of them as fundamentals and others as anti-fundamentals. 

"^^ This naturally leads to scalar states in the 3 3 = 1 (g) 8 color-adjoint representations. This mechanism is realized 

^^ in technicolor and top-color, chiral-color, and vector-like confinement [7]. In the presence of extra dimensions, color 

,— I octet scalars emerge as low-lying Kaluza-Klein modes of the bulk gluon field [S] . 

O, At the LHC sgluon pairs will be copiously produced just through their couplings to gluons. In addition, for large 

I masses, the model dependent single sgluon production channel might be competitive [1, 2^ [5]. Available studies 

^^ include the single and pairwise production in the context of supersymmetric scenarios [U |2l [12] , GUTs [H] , extra 

pH dimensions [51 [T^] , as well as more model-independent approaches [131 [H] . Color-octet vector bosons have also been 

i_i considered [15] . Distinctive decay patterns appear through couplings to pairs of SM particles or new heavy colored 

states. The most generic signature is the decay to two quark or gluon jets, pp — )■ GG* — >■ 4 jets, possibly including 

bottom jets [16]. Subjet techniques [17] or suitable cuts on jet pair invariant masses [HI [19] have been proposed to 

handle the overwhelming QCD background. In supersymmetric models with Dirac gluinos the constraints on squark 

mixing are so weak that an essentially unconstrained squark mass matrix will lead to sgluon decays to single top (anti-) 

^r^ quarks plus a light jet, G ^ tq,iq X- Finally, highly isotropic multi-jet signatures pp — > GG* — 7> tttt — 7> 8j + 2£ + $rp 

\^ are likely for sufficiently heavy sgluons [3 [TB] . Complementary rare decays [2D] or long-lived bound states [31] are 

_^ other potential discovery modes for novel color-adjoint scalars. 

^— ^ In this paper we present a complete next-to-leading order QCD calculation of sgluon pair production at the LHC. 

^ We examine the features and quantitative impact of the QCD quantum effects on the production rates and sgluon 

. . distributions. Our results are implemented in the Madgraph framework through the dedicated MadGolem package 
for the production of new particles to next-to-leading order [22] . This tool automatically computes next-to-leading 
order QCD corrections for any heavy particle production process and will be publicly available after the current 

j_^ testing phase. The NLO sgluon distributions we compare to the matched [23l [21] results for the combined process 

C3 pp ^ GG* +]eiii. 

Using the renormalizable supersymmetric realization the gluonic QCD corrections to sgluon pair production are 
obviously well defined. Additional supersymmetric QCD corrections are suppressed by the squark and gluino masses 
and thus negligible [251 I26j . Because we are only interested in sgluon pair production we can decouple all supersym- 
metric partners except for the sgluon, retaining all benefits of a renormalizable theory. This theory can as well be 
interpreted as the relevant QCD part of an effective strongly interacting theory. 

II. SGLUON PAIR PRODUCTION TO NLO 

To compute the complete NLO corrections for sgluon pair production at the LHC we minimally extend the Standard 
Model by one additional color octet, weak singlet, electrically neutral, and complex scalar field G. The sgluon couples 
to the Standard Model through the covariant derivative, D^ G^ = d^G^ + Qs f^^^ G^ A^, where A*^ denotes the 
gluon field, gs the strong coupling constant, and f^^'~" the adjoint SU{i) generators. The Feynman diagrams for the 
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Figure 1: Leading order Feynman diagrams for sgluon pair production via quark-antiquark annihilation and gluon fusion. 

two partonic LHC production mechanisms 

qq^GG* and gg ^ GG* (1) 

are shown in Figure IT] The sghion couphng to gluons reads 

C D D^G* D^G - mlGG* 

3 _ 5^ fABC [G^*(aM gB) - (^A. G^*)G^ ] Af^ + gl [/-^^^ /^^^ + /^^^ f^'^^] G^* G'^A^ A^^ . (2) 

Incidentahy, we notice the absence of direct sgluon couphngs to matter. In supersymmetry they only arise as effective 
dimension-5 operators induced by the one-loop squark and gluino loops with a non-trivial scaling for individual heavy 
masses [1^ . As long as these couplings are small — which is true if they are loop- induced — the sgluon mass range is 
not constrained by stringent bounds from dijet resonance searches. As a consequence, sgluons can be relatively light. 
Conversely, for 0(1) sgluon-quark-quark couplings sgluon masses below m^ — 0{2 TeV) are ruled out by the LHC 
experiments. 

We implement the couplings shown in Eq.S into the Madgraph framework [57]. MadGolem then generates 
all tree- level diagrams and the corresponding helicity amplitudes, making use of Madgraph and Helas [55]. All 
one-loop amplitudes and the corresponding helicity amplitudes we generate with a modified version of Qgraf [29] and 
Golem [HHl 131] • The model specific ultraviolet counter terms are part of the model implementation. The subtraction 
of infrared and (if applicable) on-shell divergences is automized [5S1 132] • 

Throughout our analysis we use the CTEQ6L1 and CTEQ6M parton densities [33] with consistent values of as. 
For the central renormalization and factorization scales we choose the average final state mass f/* = ^r.f — itig, 
which has been shown to lead to stable perturbative results [55], f5SJ . The LHC center of mass energy is \/~S = 8 TeV. 
Unless stated otherwise, we set the sgluon mass to mg = 500 GeV. 

Technically, there is a MadGraph4 |27j issue with the color structure of the quartic gluon-gluon-sgluon-sgluon 
coupling shown in Eq.S. We therefore generate the required structure f^'^^ jBDE^jAde jBCE through an auxiliary 
massive, color-adjoint vector boson V^ with an appropriate coupling to a gluon and a sgluon, namely 

GA*^B yC^. . g^ ^^ J ABC qA j^B yCi. . „^^ ^^ J ABC (3) 

The quartic gluon-sgluon interaction is then simply given by the decoupling limit n^ ^ s. For Madgraph5 this 
technical complication is not necessary any longer. 

Production rates to Next-to-Leading order 

As a first step we present the results for the total NLO cross section for sgluon pair production. Later in this 
section we focus on more specific aspects of the real and virtual corrections. Unless stated otherwise, we assume 
TOG = 500 GeV and \/5 = 8 TeV. 

The size of the QCD quantum effects we describe in terms of the consistent factor K = a^^^ ja^'^ . From the 
production of supersymmetric particles it is well known that for LHC energies of 8 TeV and particle masses in the 
500 GeV to 1 TeV mass range this correction factor can become unexpectedly large. This is not a sign of poor 
perturbative behavior but an artifact of the LO CTEQ parton densities "331 |34j . Correspondingly, Table IT] typically 
shows K > 1.5 for this collider energy while the 14 TeV scenario has smaller, yet sizeable, corrections. 

We then provide a comprehensive analysis of the LO and NLO cross sections aipp — J> GG*^ as a function of the 
sgluon mass in Fig. [2| In the left panel we show the LO and NLO cross sections with the envelope of the NLO 
factorization and renormalization scale variation in the range /i°/2 < [iR,F < 2/^". The effects of an independent as 
well as diagonal variation of the factorization and renormalization scales we show in Fig. [Sj and discuss below. As 
alluded to, the LO parton densities drive the LO cross sections to unexpectedly small values which makes the NLO 
corrections appear larger than '^ 100% for sgluon masses in the TeV range. 
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Table I: Total pp — > GG* cross sections and corresponding AT-factors for different sgluon masses and LHC energies. 

In the right panel of Fig.[2]we separate the contributions stemming from the different partonic sub-channels: qq, gg 
and the crossed purely NLO gq/gq initial state. The NLO corrections steadily increase for increasing sgluon masses. 
In part, we can trace back this behavior to threshold effects which we will further discuss below. The right panel 
of Fig. [2] also quantifies the dominance of the gluon-fusion mechanism gg — >■ GG* . The reason is twofold: first, the 
color charges in the four-octet interaction are larger than the triplet-octet combination that drives the qq channel. 
Second, the gg channel benefits from particular kinematic features of the parton level process. While qq — >■ GG* 
to LO proceeds merely through the (derivative) gGG coupling, the gg mechanism also receives a contribution from 
the (contact) quartic interaction. The first case corresponds to a p-wave and implies that the total (partonic) cross 
section scales as aqq ~ /3^ at threshold, where /3 = -y/l — Aitlq/s denotes the sgluon velocity in the center-of-mass 
frame. In the gluon fusion case the s-wave component from the quartic interaction translates into a linear dependence 
(j{gg) ~ /3. The latter dominates in the vicinity of the threshold which, moreover, corresponds to the low- a; region 
where the gluon parton densities becomes large. 

Following these arguments we can compare sgluon pair production to stop pair production (or squark pair production 
with decoupled gluinos) [26 . The differences at leading order can be traced to the relative strength of the color 
interactions arising from the fundamental vs adjoint final-state scalars. The ratios of the expected sgluon versus stop 
pair production rates can be roughly inferred from their parton level cross sections. We can compute these ratios 
directly from the corresponding analytical expressions [U [5] , 



a{qq -^ it*) 
a{qq-^GG*) 



= 1/6 



aigg^ii*) 
a{gg~^GG*) 



1/20 



(4) 



These estimates nicely agree with the NLO results for stop pair production available from Prospino [26], which give 
(j{pp — >■ it*) ~ 3 pb for stop masses of m^ ~ 350 GeV, i.e. a factor of 0(20) below the sgluon results in Table p] and 
Fig. [2] The NLO effects to squark pair production are comparatively mild. In contrast, gluino pair production as an 
example of a color-octet interaction also shows large K factors and a very pronounced dependence on the mass of the 
produced heavy particles. 
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Figure 2: LO and NLO cross sections a-{pp — >■ GG*) as a function of the sgluon mass. The band corresponds to a scale variation 
/x°/2 < fiB„F < 2^^\ In the right panels we explicitly separate the contributions from the different partonic sub-channels, qq, 
gg and also gq. 




Figure 3; Renormalization and factorization scale dependence. The plot traces the scale dependence following a contour in the 
HF'tiR plane in the range /i = (0.1 — 10) x /x° with ^^ — mc- The sgluon mass we fix to mo = 500 GeV. 

Scale dependence 

Aside from the often positive corrections to the production rate the main motivation for the computation of higher 
order corrections is the reduced theoretical uncertainty. While we cannot derive the uncertainties arising from unknown 
higher orders in QCD from first principles, we can attempt to track them in the dependence of unphysical parameters 
introduced by the perturbative approach. An example for such a parameter are the factorization and renormalization 
scales which we introduce when we remove collinear and ultraviolet divergences order by order in perturbation theory. 
In the limit of infinitely many terms in the power series in a^ these scale dependences have to vanish, so unless there 
is a systematic shift from one perturbative order to the next the scale dependence should cover the asymptotic cross 
section values. While for Drell-Yan-type processes we know that this argument fails, purely color mediated processes 
with a^'~' ex a^ tend to give a reasonable error estimate this way [551 [IS] ■ Vice versa, we can at least firmly state that 
the scale variation gives a minimum uncertainty simply because we have the freedom to choose the two scales within 
a reasonable energy range. 

In Fig. [3] we show the scale dependence for the pp — >■ GG* production rate, independently changing the renormal- 
ization {fiR) and the factorization {(J-f) scales. We illustrate the path in two dimensions in the little square in the 
first panel. The individual scale variation is chosen as fx^^^lO < /i < lO/i'-'^', where /i*^°-' stands for our central value 
choice lip — iiR = fi^^' = TUG = 500 GeV. The stabilization of the mentioned scale dependence manifests itself as 
a smoother ct^^*-* slope, with varies typically around Aa^^'-' / a^^'^ ^ 0(30%), while for leading order this variation 
can be as large as 0(80%). 

We also see that the maximum rate at small scale values which is often interpreted as a sign of scale stabilization 
is an artifact of identifying the two scales. An independent variation gives the largest rates at small values of the 
renormalization scale combined with larger values of the factorization scale — even though from a resummation point 
of view it is not clear how such a scale choice would be interpreted [Hj . 



Real emission 



Real emission corrections to sgluon pair production arise at order a^ and originate from the three-particle final state 
contributions, wherein one extra gluon accompanies the produced sgluon pair. Wc show sample Feynman diagrams 
in Fig. HI Following the standard procedure we subtract infrared divergences from the emitted gluon in the soft 
and/or collinear regimes using the massive Catani-Seymour dipoles |351 136] . In addition to the SM dipoles available 
in the MadDipole [37] package, MadGolem includes for example the sgluon dipoles to cope with our novel infrared 
divergent structure. Such divergences appear when the external sgluons radiate soft gluons and require new final-final 
and final-initial dipoles. The sgluon can also be a heavy spectator parton, but for this case we can simply use the SM 
initial-final dipoles. This is because for the dipole function the spectator carries information about the mass of the 
colored particle, but not about its spin. 

In Appendix IA] we give the new sgluon dipoles including the FKS-style phase space cutoff < a < 1 ,38^ . The 
numerical implementation is publicly available upon request. Among several numerical improvements, the parameter 
a gives us an easy handle to check our implementation. For a wide range a = 10° — 10~® we find stable cross section 
for the combination of the real emission diagrams with the integrated dipoles. As a default value in MadGolem we 
use a = lO^''. 





Figure 4: Sample Feynman diagrams for real emission corrections to sgluon pair production via quark-antiquark annihilation 
(upper) and gluon fusion (lower). 

Virtual corrections 

Virtual corrections to sgluon pair production appear as order a^ contributions from virtual gluons coupling to 
quarks and sgluons. All divergences we regularize in n = 4 — 2e dimensions. The infrared poles are cancelled after 
we include the integrated Catani-Seymour dipoles |351 I36j and take into account the collinear higher order correction 
consistently included in the definition of the parton densities. The ultraviolet divergences are absorbed in the physical 
renormalization of the strong coupling constant and the sgluon mass. As described in Appendix \m we use the MS 
scheme with decoupled heavy colored states [39, for the strong coupling and the on-shell scheme for the mass. For an 
internal check we use an independent implementation of our sgluon model in FeynArts; all MadGolem results can 
then be numerically compared to the output from FeynArts, FormCalc and LoopTools [40 . 

Starting with the dominant gluon fusion channel, in the left panel of Fig. [6] we examine different contributions to 
the real and virtual NLO corrections to the hadronic process pp — > GG* as a function of the sgluon mass. Leaving 
aside gauge invariance issues while applying a numerical test we separately show different one-loop pieces normalized 
to the LO rate, Aa^^'^ / cr^'^ . In addition, we distinguish the partonic subprocesses qq and gg. The crossed channel 
qg does not develop any virtual corrections but is required for the complete cancellation of the collinear divergence. 

The real emission together with the virtual box diagrams contributes the bulk of the NLO quantum effects. Both 
feature a characteristic growing trend with increasing sgluon mass. For intermediate sgluon masses, real emission 
gives rise to corrections in the ball-park of 20 — 60%, but it may eventually reach up to 100% for TeV-scale sgluons. 
Note that this behavior cannot be interpreted as a break-down of perturbation theory because QCD corrections to 
our supersymmetric setup are fundamentally well defined. Box-like topologies, i.e. one-loop corrections to the ggGG 
vertex including diagrams shown as the 3th and 4th diagrams in the lower row of Fig. [5J amount to roughly 40% 
and exhibit a slightly milder dependence on ttiq. Both, the size of these contributions and their increase with mc 
we can attribute to the peculiar threshold behavior of the NLO corrections. In particular, long-range gluon exchange 
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Figure 5: Sample Feynman diagrams for virtual corrections to sgluon pair production via quark-antiquark annihilation (upper) 
and gluon fusion (lower). 
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Figure 6: Relative size Aa'^^^/a 



NLO /„LO 



a )/a of the real emission and virtual corrections to cr{pp — >■ GG") as a 



function of the sgluon mass ma- We separate the partonic gg (left) and qq (right) initial states. The contribution from the 
self-energies is negligible and not explicitly shown. 



between slowly moving heavy final-state particles /3 — > develops a Coulomb singularity a^^'~' ^ T^cts/ P which cancels 
the linear dependence from the tree-level contribution a^'^ ^ (3 and leads to a finite NLO rate but a divergent K 
factor [26j. Gluon radiation off the initial state carries positive, and potentially large, logarithmic pieces which supply 
an additional source of enhancement - and that can eventually be resummed [d\ . 

For the subleading qq initial state we find a sizable contribution from the giuon-sgluon-sgluon gGG vertex correc- 
tions. They are fairly independent of variations of the sgluon mass. The remaining one-loop topologies only contribute 
at the percent level of less. 



III. DISTRIBUTIONS: NLO VERSUS MULTI-JET MERGING 



Predictions based on (fixed-order) NLO cross sections entail significant improvements of the central values and the 
theory uncertainties, as we have just shown for the specific case pp — ?> GG* . Before these results can be integrated 
in experimental analyses we need to confirm that this quantitative picture also holds for the main distributions. 
Previous work in the literature shows that the transverse momentum and rapidity distributions of pair produced 
heavy particles are relatively stable with respect to higher-order corrections 22i i25i . For the production of such 
heavy particles a parton shower should deliver a good description of jet radiation patterns, because its underlying 
coUinear approximation applies for a wide range of transverse momenta relative to the sgluon masses (TJ 1411 I42j . 
Normalizing the event numbers generated by standard Monte Carlo tools to the NLO cross section should therefore 
be an appropriate strategy. 

To quantitatively assess such statements we compare the fixed order NLO parton-level distributions for the produc- 
tion process pp -^ GG* , as obtained from MadGolem, with a multi-jet merging calculation. For the latter we employ 
the Mlm [53] scheme and generate events using MadGraph 4.5 [32] interfaced with Pythia [13]. The entire descrip- 
tion of our sgluon is supplied by the MadGolem model file. In Fig. [7] we display the resulting transverse momentum 
and rapidity distributions for one outgoing sgluon for the NLO calculations as well as for jet merging including up 
to two hard jets. Only including one hard jet in the merging prescription would not change the results within their 
numerical precision. Those two distributions are normalized to unity. For the NLO results we separately show the 
LO, real emission, and virtual gluon contributions, defined in terms of Catani-Seymour dipoles with a = 10"'^. 

First, we see that all different ingredients of the NLO distributions have essentially the same shapes. Large effects 
on the total rate from collinear radiation or Coulomb singularities only have a negligible effect on the distributions of 
the heavy states. In addition, the normalized fixed-order and merged distributions agree very well. Small differences 
like the slightly harder px profile of the merged prediction are accounted for by the extra recoil jets from the parton 
shower. Similarly, such a second jet from initial state radiation can balance the first emission and lead to more central 
sgluons in the detector. 
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Figure 7: Sgluon transverse momentum and rapidity distributions at parton level. We assume ma ~ 500 GeV and \^ = 8 TeV. 
For the NLO curves we separately display the LO, virtual, and real contributions (a = 10~^). In addition, we show the 
corresponding distributions based on multi-jet merging in the MLM scheme [23] with up to two hard radiation jets. The NLO 
and merged results are normalized to unity while the different contributions to the NLO rates are shown to scale. 

IV. SUMMARY 

We report on the first complete calculation of sgluon pair production to next-to-leading order. The pairwise 
production of scalar color-adjoints we define in terms of an extended supersymmetric model; however, after decoupling 
squarks and gluinos our results can be considered reasonably model-independent. The sgluons have tree level QCD 
couplings to gluons and do not couple to matter. We find 

1. potentially large production rates, driven by the gluon fusion subprocess. Typical numbers range around a '^ 
0(1 pb) for sgluon masses of rriG — 500 GeV and \/S — 8 TeV. Such small masses are not ruled out by current 
experimental constraints. 

2. sizable NLO quantum effects, traced back to real gluon radiation and a certain subsets of vertex and box virtual 
corrections. Their relative size increases with the sgluon mass, mainly due to threshold effects. 

3. substantially reduced theoretical uncertainties. The scale dependence which is dominated by the renormalization 
scale choice and which may reach 0(80%) at leading order, is reduced by a factor 1/2 — 1/4. 

4. NLO sgluon distributions which agree very well with complementary results from multi-jet merging. Applying 
the NLO rate normalization to sgluons-|-jets event generation should give very consistent predictions for the 
LHC. 

Besides its phenomenological impact our study illustrates the performance of the (soon-to-be-public) MadGolem 
package. The genuine dipole and counter term structures which cope with infrared and ultraviolet divergences in the 
presence of the sgluon field have been implemented and can be accessed automatically. In this sense the present study 
qualifies as a non-trivial example of a fully automized NLO calculation for the production of heavy particles beyond 
the Standard Model. 
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Appendix A: Sgluon dipoles 



Sgluons are color octets with spin zero, so their dipoles are identical to supersymmetric scalar quarks with modified 
color factors Cp — >■ Ca- To remove the related infrared divergences we implement the (un)integrated dipoles presented 
in Appendix C of Ref. |36) with this replacement. In addition, we introduce a variable size of the subtraction phase 
space, < a < 1, as pioneered in the FKS subtraction scheme [55] • Values a < 1 limit the phase space region over 
which we subtract finite dipole contributions around the soft-coUinear pole. Our notation closely follows Ref. [36] . 

From Eq.(C.l) of Ref. [5S] we obtain the sgluon dipole function (VgCk) for the final-final case. The corresponding 
integrated dipole is decomposed into an eikonal part including the soft integrals and the remaining hard collinear 
integrals. 
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PniPj,Pk) = 



1 + Vij^k 



Pi 



1 + u,j- fe + 2^2 / (1 - p] - pI) 

The kinematics of the splitting is described by 

PtPk 



[n = j, k) . 



= 1- 



Uij.k 



PiPj 



PiPk + PjPk 

where the upper limit in the jjij^k phase space integration is 

2pk (1 - Pk) 



PiPj + PiPk + PjPk 



2/+ = l- 



^-Pf 



PI 



Pi 



(A3) 



(A4) 



(A5) 



To include the phase space parameter a in the unintegrated dipole we simply add the step function Q{a — ygG,kly+)- 
It ensures that only for HgCk < ot y+ instead of the entire range ygG,k < y+ the dipole is subtracted from the hard 
matrix element. For the integrated dipole part we modify Eq.( Al ) by subtracting the finite phase space contributions 
which a choice of a ^ 1 removes. 



IgG.k (a) 



= /, 



IgG,k - ^IgG.k [o) 

27r 
as 



gG,k 



[dPg iPgG,Pk)] 6 ( -^^ - a 



{VgG,k) 

2pgPG 



(A6) 



For the (by definition) finite contribution AlgQ^ we set e = 0. The eikonal part of the kernel 2/[l — zq {1 — HgCk)] is 
the same for the {VgQ^k) and {VgG,k), so we use the SM result for AI°^^{a) as provided in Eq.(A.9) of Ref. 44 . The 
collinear part is different, giving a correction to Eq.( Al I of the form 



- AIgG,k (a) = -Ca 



2AI 




(1-M/c) 



1 



Pg 



^(l-a)+loga 



Pk 



(A7) 
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Figure 8: a dependence of the final-final (left) and final-initial (right) sgluon dipoles for the sub-process gg — )■ GG* g. 



As expected, this result becomes trivial for a = \ and diverges in the limit a — > 0. 

For the final-initial dipole function we start from Eq.(C.3) of Ref. [36 which gives rise to the regularized integrated 
dipole function 



Ila (x) = Ca [j^g (^)+ + 5 (f - a:) ( j;^^ + j;^^^) 



(A8) 



The three contributions to I^q are 



J^G (^)+ = 



'-2-21og(l- 



x + ^il) 



l-x 



1-x 



log {2 + fil- 



Ta;S _ 1 r,^ -2e { ^ ,^ 



e 



J;g' -2C2 - 2Li2 



f 



1 + Mg 



C2 + 2 j - - log (f + /4) 
2Li2(-M^)-^log2(f + M^) 



-4-C2 



(A9) 



Again, we introduce an a parameter into the unintegrated phase space integration, limiting the application of the 
dipole subtraction to the region f — Xgc.a < o- The kinematical variable Xij,a is given by 



^tj,a 



PaPi + PaPj - PiPj + \ {mfj - mf - mf) 

PaPi + PaPj 



(AlO) 



The additional contribution to the integrated dipole I^q shown in Eq. ( A8 ) is 



AI^G (a) ^ C. ^(\--) (-2 + 21og (1 + _; 



(All) 
Again, this result becomes trivial for a — 1 and diverges in the limit of a — > now when performing the integral over 



The numerical effects of a variable a parameter on the subprocess gg — >■ GG* g we show in Fig. [8| While the 
individual unintegrated and integrated dipole contributions diverge logarithmically with small a the sum of them is 
numerically stable over eight orders of magnitude. This kind of test is sensitive to many aspects of our calculation, 
namely the proper coverage of all divergences, the relative normalization of the two and three particle phase space, etc. 
In Fig. ^ we see that a default value of a = fO~^ gives roughly equal contributions from unintegrated and integrated 
dipoles, avoiding large numerical cancellations for the final-initial dipole. For the final-final dipole we would have to 
go to smaller values of a which make the final-initial case harder, so we use a — IQ^^ throughout. 



Appendix B: Renormalization 



The ultraviolet counter terms we include automatically via the leading-order Qgraf output. At present, Mad- 
GOLEM fully supports the calculation of NLO QCD corrections for the Standard Model, the MSSM, and several 
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G 

Figure 9: Feynman diagrams for the sgluon field renormalization (left) and sgluon-mediated gluon field renormalization (right). 

other extensions of the Standard Model, including sgluons. For dimensional regularization we employ the standard 't 
Hooft-Veltman scheme with n = A — 2e dimensions. The renormalization constants we define through the additive or 
multiplicative relations between the bare and the renormalized fields, 

*(")-> Z^^* m'-^^ ^ m^ + Sm^ gf^^gs + 6g, (with * = q,A,G). (Bl) 

These field, mass and coupling renormalization constants we conventionally phrase in terms of two-point functions 
which we supply in a separate library. Given a generic Lagrangian /^(^P,™^, g^) with a QCD interaction this consis- 
tently gives a counter term Lagrangian of the form (5£(^, mq,, g^, S'^ , Sm^j,, Sgs). 

First of all, a new sgluon field modifies the strong beta function. If we start with the quantum corrections to the 
quark-quark-gluon vertex in terms of the strong coupling Zg_, , the gluon field renormalization Z3 , and the quark field 

1/2 
renormalization Z2 this translates into a combined Zi = Zg^ ^2 Z^ . Each of these renormalization constants we 

expand as Zi = 1 + Si + 0{a'^), with MS counter terms Si. The strong coupling constant renormalization at one loop 

we can thus write as 

Sgs ^ Si- S2- -S3 



2 



rSM '^s 



with Si =St''' = - -^ (Ca + Cf) Ae 

47r 



S^^Sl'^^-^CpA, 



in 



Ss ^Sl'^+Si = ^ (^^CA-ufC^Tn^ _ _|£_ C^ A, . (B2) 

The shifted pole in the MS prescription is A^ — \/e — 7s + log(47r) and the active number of fermions is Uf — 6. 
Because there are no direct couplings between sgluons and matter fields S2 keeps its SM value. For the same reason, 
sgluon-mediated corrections to the quark-quark-gluon vertex are absent at one loop, so Si does not change. Only the 
gluon self energy is modified by the triple and quartic gluon/sgluon interactions, as displayed in Fig. ^ 

Combining all of the above contributions and decoupling the heavy {H) colored degrees of freedom — in our case 
the top and the sgluon — gives us the final expression for Sgs in terms of the measured as values. We implement 
this subtraction in the zero-momentum scheme [251 I39j . It leaves the renormalization group running of a^ merely 
determined by the light (L) degrees of freedom. The renormalization constant finally reads 

_ asP^+P^. asfl.^ml I ml 



A, - -^ - log ^ + - log 



2 ' o ^"6 2 



47r 2 " 47r V3 ° /i|j 2 ' ^^ 

/3o = P^ + P^ ^(^-^Ca- {uf - 1) Cf Tn\ - fcp T^ + ^ Ca) ■ (B3) 

In a second step we need to compute the QCD renormalization constants in the sgluon sector. The sgluon two-point 
function receives 0{as) corrections due to virtual gluon interchange, as shown in Fig.^ The corresponding ultraviolet 
divergences we absorb into the sgluon mass mo and field-strength Zq. As renormalization condition we choose the 
on-shell scheme 

^et'imD^O => 5Zg = - 5Re E'(m^) 
3ffeE(TO^)=0 => (5tog= + 3fie E(m^) , (B4) 

where Jfte Eg denotes the (real part of the) renormalized sgluon self-energy, 

taiq^) = Sg(?') + (9' - ml) SZg - Sml , (B5) 
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^JUUUUx;' „ -igsf-^^^ [Sgs + ^ {SZa + SZa- + SZa)] [G*^(9^ G^) - (9^ G*-^)G^] A^ 

AC \ G" 



G' .-' 




i gl (/^^^ /''''"' + f^^'^ f^^^) [2 5ff, + SZa + SZa] G*^ G^ A;^ A^ 



p^ (5Zg — Sttiq — rriQ SZq 



Table II: Counter term Feynman rules for the sgluon-mediated interactions. 

and S'((7^) = cP /dq^ ^(9^) the corresponding derivative with respect to the momentum squared. The analytic form 
of all renormalization constants we reduce down to one and two-point scalar loop integrals _45] ■ The sgluon mass and 
field strength renormalization then reads 



5Zg ^ i^^A [Bo(r7ig,TOG,0) + m%B'f^{m%,m%,0)] 



Sma = -Ca 



"^G + ^^o("4] 



(B6) 



As expected, these expressions are identical to the squark case, modulo a factor Ca/Cf that reflects the different 
SUc{3) representations. 

Finally, in Table [IT] we quote the analytical expressions for the relevant ultraviolet counter terms SC as a function 
of the field, mass, and strong coupling renormalization constants derived in this Appendix. 
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